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C**- ! Using gauge invariant hadronic two- and three- density correlators we extract informa- 
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tion on the spatial distributions of quarks in hadrons, and on hadron shape and multipole 



>• ■ moments within quenched lattice QCD. Combined with the calculation of N to A transi- 

^ ■ tion amplitudes the issue of nucleon deformation can be addressed. 

O _ 

CN ' 1. Introduction 



Correlation functions calculated in lattice QCD can be connected to basic hadronic 
features. Two- and three- density correlators for mesons and baryons reduce in the non- 
relativistic limit to the square of the wave function and therefore provide detailed informa- 
tion on hadron structure such as quark spatial distributions, hadronic shapes and charge 
radii. In addition lattice results can be used to test predictions in various models. Direct 
connection to the issue of nucleon deformation, currently under experimental investiga- 
tion j 1 112 j . is made by calculating the 7*iV to A transition form factors in quenched and 
unquenched lattice QCD. The aim is to obtain an accurate determination of the ratio of 
the electric quadrupole to the magnetic dipole amplitudes, Rem(q 2 ), as a function of the 
momentum transfer q. A non-zero Rem arises due to deformation in the nucleon and/or 
A and it is attributed to different mechanisms in the various models: in quark models 
the deformation arises due to the colour-magnetic tensor force whereas in 'cloudy' baryon 
models due to the meson exchange currents. 

Let us first consider the equal-time correlators [3], 

C*(T,t)=jd 3 r'{H\ft(r\t)p*(r' + r,t)\H) , (1) 

with /3p(r, t) given by the normal order product : u(r, t)Tu(r, t) :. For T = 70 and T = 1 we 
obtain, in the nonrelativistic limit, the charge and matter density respectively. For T = 75 
we obtain the pseudoscalar density which we will compare to bag model predictions. The 
matrix elements of two-density correlators are shown schematically in Fig. ^ The advan- 
tage of using density correlators is that they are gauge-invariant unlike Bethe-Salpeter 
amplitudes. For baryon wave functions three-density correlators are needed. Here we will 
only show results for the one-particle density obtained after integrating the wave function 
over one relative distance. The diagram involved is shown in the lower part of Fig. ^ 

Since the spectroscopic quadrupole moment of the nucleon vanishes in order to deter- 
mine any possible deformation we evaluate the N to A electromagnetic transition form 

*Talk presented by C. Alexandre-it. 
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Figure 2. Time evolution to ground state for 
the rho meson at k = 0.153. 



Figure 1. Equal-time correlator for a 
meson (upper) and abaryon (lower). 

factors. We calculate the three-point correlation functions |JJ (G A ^ N {t 2) t; p p; T)) and 
{G^' lA (t 2 , t; p p; T)) and consider the following ratio |3] 

■ (G A ^ N (t 2 , t; p p; r)> (G^ A (t 2 , t- -p, -p T^ 1 ' 1 



t 2 -t > l,t > 1 



(-g l3 G AA (t 2 , p r 4 )> (G™(t 2 , -p; r 4 )> 
E N (p) + M N y/2[E' A (p')+M N 



2E N {p) 



2E> A (p>) 



1/2 _ 



(2) 



from which the transition form factors can be extracted. The index a is the Lorentz 
index for a spin- | field, j^(x), is the lattice conserved electromagnetic current inserted 
at time t, E^(p) and E' A (p') are the energies of the nucleon and the A respectively and 
t 2 , the time location of the sink, is varied to identify the ground state. We use projection 
matrices Tj = 1/2^°) and r = 1/2 (o o) as m re f- 0- 
2. Results 

We have analysed 220 quenched configurations at f3 = 6.0 for a lattice of size 16 3 x 32 at 
k= 0.15, 0.153, 0.154, 0.155 which give ratio ofpion to rho mass of 0.88, 0.84, 0.78, 0.70 
respectively. Using the standard definition of the naive quark mass, 2m q a = (1/k — 1/k c ), 
where k c is the value of k at which the pion becomes massless, we find m q ~ 300, 170, 130 
and 85 Mev respectively. We used the value of the string tension to set the physical scale 
obtaining for the inverse lattice spacing a -1 ~ 1.94 GeV. The density insertions are taken 
at t = T/4 = 8a. To check that this time interval is long enough to extract the physical 
hadronic ground state we computed the density-density correlators for varying values of 
the insertion time t. Fig. |21 shows the results for the rho correlators which converge when 
t > 6a. Moreover, we note a clear deformation (ratio C^ Q (x, 0, 0)/C^ o (0, 0, z) ^ 1) which 
remains approximately the same as early as t = 3a even though the z- and x-profiles 
change appreciably showing that the deformation is a robust, physical property of the rho 
meson in its ground state as well as its low-lying excited states. 

The dependence of hadron wave functions on the quark mass is displayed in Fig. El 
where we plot the density-density correlators at various k values for the pion, the rho, the 
nucleon and the A + normalized to unity. The size of the nucleon and A wave functions 
do not change for naive quark mass smaller than 130 MeV. The rho wave function shows 
the largest variation with the quark mass. In FigEJwe compare the matter and charge 
densities at k — 0.153. We find that the matter density falls off more rapidly than the 
charge density and it has the same shape for all four hadrons. 
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Figure 3. Density-density correlators as 
function of the quark mass. 
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a Figure 4. Comparison of charge (a) and 
matter (b) densities at k = 0.153. 
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Figure 5. The pseudoscalar density for 
the pion at k — 0.153,0.154 and 0.155 
with fits to r 2 (a + br 2 + cr 4 ) exp(— mr). 



Figure 6. Pseudoscalar density for the pion 
(upper), the rho (middle) and the nucleon 
(lower) at k = 0.153. Fits as in Fig. 



The pseudoscalar density for the pion is shown in Fig. El for m^jm p = 0.84, 0.78 and 
0.70. In Fig. |U]it is compared with that for the rho and the nucleon at m n /m p = 0.84. The 
bag model predicts that the integral / d 3 rC^ 5 (r) is zero 0. A reasonable fit is obtained 
by an exponential times a polynomial ansatz. The long tail of the data and the integral 
of the fitted ansatz both favor a non-zero integral. However a careful extrapolation to 
the continuum limit (using large volumes) is required for this quantity, especially with 
Wilson fermions as used here. 

The rho deformation seen in Fig. El can be made more quantitative by analysing the 
density-density correlator into a dominant L = state and a suppressed L = 2 [7] : 



<Pi(0)|p 70 (r)/5 70 ( 0)|p i (0) >= 



(3 X 2 _ r 2\ 

>l( r ) + 3 3r2 Mr) 



(3) 



where \pj(0) > is a zero momentum state with polarization j. 
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Figure 7. Decomposition of C 70 for the rho into 



As shown in Fig. a good descrip- 
tion of the rho correlator, is obtained 
by taking <p±(r) = Aexp(— m\r) 
and 02( r ) = Br 2 exp(— rri2r). Us- 
ing the values extracted from the fits 
and neglecting B 2 terms we find for 
the deformation 5 = (3/4) ((3z 2 — 
r 2 )/(r 2 )) ~ 0.01 with an error of 
about 80% which mainly arises from 
the poor determination of the coeffi- 
cient of the L = 2 state. A direct de- 
termination of the quadrupole mo- 
ment from the rho correlator yields 
5 = 0.03 ± 0.01 in reasonable agree- 
ment with the value obtained from 
angular momentum part L = and L = 2. the angular decomposition. 

First lattice results for the ratio of the electric quadrupole to magnetic dipole ampli- 
tudes, Rem, for both quenched and for two dynamical Wilson fermions are obtained at 
the smallest allowed lattice momentum, q = (2ir/L s , 0, 0), where L s is the spatial length 
of the lattice. At q 2 = 0.54 GeV 2 for the unquenched theory using the SESAM configu- 
rations jH] we find in the chiral limit Rem ~ (—3.7 ± 0.5)%. For q 2 = 0.15 GeV 2 in the 
quenched theory on a lattice of size 32 3 x 64 at (3 = 6.0 we find Rem ~ (—3.0 ± 0.3)%. 
We used the nucleon mass to set the physical scale. Finite lattice spacing effects on the 
Rem ratio can be significant and are currently under investigation. The finite volume 
dependence on Rem is also under study. 
3. Conclusions 

We have presented a gauge-invariant determination of hadron profiles in the quenched 
approximation. We have found that the rho wave function shows the strongest dependence 
on the quark mass whereas the nucleon and A the weakest. We have established that the 
rho is deformed with deformation which increases as we approach the chiral limit, whereas 
the A + has no statistically significant deformation. The matter density is similar for all 
four hadrons and falls off more rapidly than the charge density. Using state-of-the-art 
lattice techniques the phenomenologically important Rem has been extracted for two q 2 
values and found to be consistent with experimental measurements. 
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